skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Search for: All records

Creators/Authors contains: "Stock, Charles"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. Abstract. Zooplankton diel vertical migration (DVM) is critical to ocean ecosystem dynamics and biogeochemical cycles, by supplying food and injecting carbon to the mesopelagic ocean (200–800 m). The deeper the zooplankton migrate, the longer the carbon is sequestered away from the atmosphere and the deeper the ecosystems they feed. Sparse observations show variations in migration depths over a wide range of temporal and spatial scales. A major challenge, however, is to understand the biological and physical mechanisms controlling this variability, which is critical to assess impacts on ecosystem and carbon dynamics. Here, we introduce a migrating zooplankton model for medium and large zooplankton that explicitly resolves diel migration trajectories and biogeochemical fluxes. This model is integrated into the MOM6-COBALTv2 ocean physical-biogeochemical model, and applied in an idealized high-resolution (9.4 km) configuration of the North Atlantic. The model skillfully reproduces observed North Atlantic migrating zooplankton biomass and DVM patterns. Evaluation of the mechanisms controlling zooplankton migration depth reveals that chlorophyll shading reduces by 60 meters zooplankton migration depth in the subpolar gyre compared with the subtropical gyre, with pronounced seasonal variations linked to the spring bloom. Fine-scale spatial effects (<100 km) linked to eddy and frontal dynamics can either offset or reinforce the large-scale effect by up to 100 meters. This could imply that for phytoplankton-rich regions and filaments, which represent a major source of exportable carbon for migrating zooplankton, their high-chlorophyll content contributes to reducing zooplankton migration depth and carbon sequestration time. 
    more » « less
  2. Abstract Phytoplankton stoichiometry modulates the interaction between carbon, nitrogen and phosphorus cycles. Environmentally driven variations in phytoplankton C:N:P can alter biogeochemical cycling compared to expectations under fixed ratios. In fact, the assumption of fixed C:N:P has been linked to Earth System Model (ESM) biases and potential misrepresentation of responses to future change. Here we integrate key elements of the Adaptive Trait Optimization Model (ATOM) for phytoplankton stoichiometry with the Carbon, Ocean Biogeochemistry and Lower Trophics (COBALT) ocean biogeochemical model. Within a series of global ocean‐ice‐ecosystem retrospective simulations, ATOM‐COBALT reproduced observations of phytoplankton N:P, and compared to static ratios, exhibited reduced phytoplankton P‐limitation, enhanced N‐fixation, and increased low‐latitude export, improving consistency with observations and highlighting the biogeochemical implications of dynamic N:P. We applied ATOM‐COBALT to explore the impacts of different physiological mechanisms hypothesized to underlie N:P variation, finding that two mechanisms together drove the observed patterns: proportionality of P‐rich ribosomes in phytoplankton cells to growth rates and reductions in P‐storage during scarcity. A third mechanism which linked temperature with phytoplankton biomass allocations to non‐ribosomal proteins, led only to relatively modest impacts because this mechanism decreased the temperature dependence of phytoplankton growth rates, compensating for changes in N:P. We find that there are quantitative response differences that associate distinctive biogeochemical footprints with each mechanism, which are most apparent in highly productive low‐latitude regions. These results suggest that variable phytoplankton N:P makes phytoplankton productivity and export resilient to environmental changes, and support further research on the physiological and environmental drivers of phytoplankton stoichiometry and biogeochemical role. 
    more » « less
  3. Abstract. As a key biogeochemical pathway in the marine nitrogen cycle, nitrification (ammonia oxidation and nitrite oxidation) converts the most reduced form of nitrogen – ammonium–ammonia (NH4+–NH3) – into the oxidized species nitrite (NO2-) and nitrate (NO3-). In the ocean, these processes are mainly performed by ammonia-oxidizing archaea (AOA) and bacteria (AOB) and nitrite-oxidizing bacteria (NOB). By transforming nitrogen speciation and providing substrates for nitrogen removal, nitrification affects microbial community structure; marine productivity (including chemoautotrophic carbon fixation); and the production of a powerful greenhouse gas, nitrous oxide (N2O). Nitrification is hypothesized to be regulated by temperature, oxygen, light, substrate concentration, substrate flux, pH and other environmental factors. Although the number of field observations from various oceanic regions has increased considerably over the last few decades, a global synthesis is lacking, and understanding how environmental factors control nitrification remains elusive. Therefore, we have compiled a database of nitrification rates and nitrifier abundance in the global ocean from published literature and unpublished datasets. This database includes 2393 and 1006 measurements of ammonia oxidation and nitrite oxidation rates and 2242 and 631 quantifications of ammonia oxidizers and nitrite oxidizers, respectively. This community effort confirms and enhances our understanding of the spatial distribution of nitrification and nitrifiers and their corresponding drivers such as the important role of substrate concentration in controlling nitrification rates and nitrifier abundance. Some conundrums are also revealed, including the inconsistent observations of light limitation and high rates of nitrite oxidation reported from anoxic waters. This database can be used to constrain the distribution of marine nitrification, to evaluate and improve biogeochemical models of nitrification, and to quantify the impact of nitrification on ecosystem functions like marine productivity and N2O production. This database additionally sets a baseline for comparison with future observations and guides future exploration (e.g., measurements in the poorly sampled regions such as the Indian Ocean and method comparison and/or standardization). The database is publicly available at the Zenodo repository: https://doi.org/10.5281/zenodo.8355912 (Tang et al., 2023). 
    more » « less
  4. Abstract. In marine ecosystems, most physiological, ecological, or physical processes are size dependent. These include metabolic rates, the uptake of carbon and other nutrients, swimming and sinking velocities, and trophic interactions, which eventually determine the stocks of commercial species, as well as biogeochemical cycles and carbon sequestration. As such, broad-scale observations of plankton size distribution are important indicators of the general functioning and state of pelagic ecosystems under anthropogenic pressures. Here, we present the first global datasets of the Pelagic Size Structure database (PSSdb), generated from plankton imaging devices. This release includes the bulk particle normalized biovolume size spectrum (NBSS) and the bulk particle size distribution (PSD), along with their related parameters (slope, intercept, and R2) measured within the epipelagic layer (0–200 m) by three imaging sensors: the Imaging FlowCytobot (IFCB), the Underwater Vision Profiler (UVP), and benchtop scanners. Collectively, these instruments effectively image organisms and detrital material in the 7–10 000 µm size range. A total of 92 472 IFCB samples, 3068 UVP profiles, and 2411 scans passed our quality control and were standardized to produce consistent instrument-specific size spectra averaged to 1° × 1° latitude and longitude and by year and month. Our instrument-specific datasets span most major ocean basins, except for the IFCB datasets we have ingested, which were exclusively collected in northern latitudes, and cover decadal time periods (2013–2022 for IFCB, 2008–2021 for UVP, and 1996–2022 for scanners), allowing for a further assessment of the pelagic size spectrum in space and time. The datasets that constitute PSSdb's first release are available at https://doi.org/10.5281/zenodo.11050013 (Dugenne et al., 2024b). In addition, future updates to these data products can be accessed at https://doi.org/10.5281/zenodo.7998799. 
    more » « less
  5. Abstract. The coastal ecosystem of the Gulf of Alaska (GOA) is especially vulnerable to the effects of ocean acidification and climate change. Detection of these long-term trends requires a good understanding of the system’s natural state. The GOA is a highly dynamic system that exhibits large inorganic carbon variability on subseasonal to interannual timescales. This variability is poorly understood due to the lack of observations in this expansive and remote region. We developed a new model setup for the GOA that couples the three-dimensional Regional Oceanic Model System (ROMS) and the Carbon, Ocean Biogeochemistry and Lower Trophic (COBALT) ecosystem model. To improve our conceptual understanding of the system, we conducted a hindcast simulation from 1980 to 2013. The model was explicitly forced with temporally and spatially varying coastal freshwater discharges from a high-resolution terrestrial hydrological model, thereby affecting salinity, alkalinity, dissolved inorganic carbon, and nutrient concentrations. This represents a substantial improvement over previous GOA modeling attempts. Here, we evaluate the model on seasonal to interannual timescales using the best available inorganic carbon observations. The model was particularly successful in reproducing observed aragonite oversaturation and undersaturation of near-bottom water in May and September, respectively. The largest deficiency in the model is its inability to adequately simulate springtime surface inorganic carbon chemistry, as it overestimates surface dissolved inorganic carbon, which translates into an underestimation of the surface aragonite saturation state at this time. We also use the model to describe the seasonal cycle and drivers of inorganic carbon parameters along the Seward Line transect in under-sampled months. Model output suggests that the majority of the near-bottom water along the Seward Line is seasonally undersaturated with respect to aragonite between June and January, as a result of upwelling and remineralization. Such an extensive period of reoccurring aragonite undersaturation may be harmful to ocean acidification-sensitive organisms. Furthermore, the influence of freshwater not only decreases the aragonite saturation state in coastal surface waters in summer and fall, but it simultaneously decreases the surface partial pressure of carbon dioxide (pCO2), thereby decoupling the aragonite saturation state from pCO2. The full seasonal cycle and geographic extent of the GOA region is under-sampled, and our model results give new and important insights for months of the year and areas that lack in situ inorganic carbon observations. 
    more » « less
  6. Abstract Numerical models of ocean biogeochemistry are relied upon to make projections about the impact of climate change on marine resources and test hypotheses regarding the drivers of past changes in climate and ecosystems. In large areas of the ocean, iron availability regulates the functioning of marine ecosystems and hence the ocean carbon cycle. Accordingly, our ability to quantify the drivers and impacts of fluctuations in ocean ecosystems and carbon cycling in space and time relies on first achieving an appropriate representation of the modern marine iron cycle in models. When the iron distributions from 13 global ocean biogeochemistry models are compared against the latest oceanic sections from the GEOTRACES program, we find that all models struggle to reproduce many aspects of the observed spatial patterns. Models that reflect the emerging evidence for multiple iron sources or subtleties of its internal cycling perform much better in capturing observed features than their simpler contemporaries, particularly in the ocean interior. We show that the substantial uncertainty in the input fluxes of iron results in a very wide range of residence times across models, which has implications for the response of ecosystems and global carbon cycling to perturbations. Given this large uncertainty, iron fertilization experiments based on any single current generation model should be interpreted with caution. Improvements to how such models represent iron scavenging and also biological cycling are needed to raise confidence in their projections of global biogeochemical change in the ocean. 
    more » « less